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We compare the numerical predictions of the A>r-factorization approach (semi- 
hard theory) for heavy quark production in high energy nucleon-nucleon and 
photon-nucleon collisions with the experimental data from Tevatron-collider and 
HERA. Predictions for heavy quark production at Tevatron, LHC and HERA are 
also presented. 
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1 Introduction 

The investigation of heavy quark production in high energy hadron collisions is an 
important method for studying the quark-gluon structure of hadrons. The description 
of hard interactions in hadron collisions within the framework of QCD is possible only 
with the help of some phenomenology, which reduces the hadron-hadron interaction 
to the parton-parton one via the formalism of the hadron structure functions. The 
cross sections of hard processes in hadron-hadron interactions can be written as the 
convolutions of sub-process QCD matrix elements squared with the parton distributions 
in the colliding hadrons. 

The most popular and technically simplest approach is the so-called QCD collinear 
approximation, or parton model (PM). In this model all particles involved are assumed to 
be on mass shell, carrying only longitudinal momenta, and the cross section is averaged 
over two transverse polarizations of the incident gluons. The virtualities q 2 of the 
initial partons are taken into account only through their structure functions. The cross 
sections of QCD subprocess are calculated usually in the next to leading order (NLO) 
HI, @, [| f|, H of as expansion. The transverse momenta of the incident partons are 
neglected in the QCD matrix elements. This is the direct analogy of the Weizsaecker- 
Williams approximation in QED. 

Another possibility to incorporate the incident parton transverse momenta is referred 
to as kr- factorization approach |], [7], |8|, 51 [10], 11], or the theory of semihard interactions 



12] , [13| , |14] , [15| , [Tq , |17| [18], [19], |2(|. Here the Feynman diagrams are calculated with 
account of the virtualities and possible polarizations of the incident partons. In the 
small x domain there are no grounds to neglect the transverse momenta of the gluons, 
qiT and q2T, i n comparison with the quark mass and transverse momenta, Pit- Moreover, 
at the very high energies and PiT the main contribution to the cross sections comes from 
the region qir ~ Pit or g 2 T ~ Pit, see [^T], EE| |2Bj for details. The QCD matrix elements 
of the sub-processes are rather complicated in such an approach. We have calculated 
them in the LO. On the other hand, the multiple emission of soft gluons is effectively 
included here. That is why it is the question which approach is more constructive. 

The A/r-factorization approach is based on so-called unintegrated parton distributions 
which, at the moment, are known with insufficient accuracy. Here we find them with the 



help of the realistic gluon distribution GRV94 ]23] which is compatible with the most 
recent data, see discussion in Ref. p5| . 

In Sect. 2 we present very shortly the formalism of the k^- factorization approach. 
The comparison of our numerical results with some experimental data on beauty pro- 
duction at Tevatron-collider together with the predictions including LHC energy are 
discussed in Sect. 3. Predictions of the /cr-factorization approach for heavy quark pho- 
toproduction for HERA are given in Sect. 4. 



Heavy quark production in the /c^-factorization 
approach 



The main contribution to the cross section of heavy quark production at small x is known 
to come from gg — > QQ subprocess (and jg — > QQ in the case of photoproduction). 
The transverse momenta of the incident gluons in the small-x region result in the k?- 
factorization approach from as In k\ gluon diffusion. This is described by the so-called 



unintegrating gluon distribution which can be written in the form [22 
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where the survival probability T g that gluon g with the transverse momentum q? remains 
untouched in the evolution up to the factorization scale \x is 
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Here a'(x, q?) denotes xg(x,q^) or xq(x,q^), and P ga '(z) are the splitting functions. 

We have to emphasize that f g (x,qT,fi) is just the quantity which enters into the 
Feynman diagrams. The distributions f g (x,qT,fi) involve two hard scales^: qr and the 
scale n of the probe. The scale [i plays a dual role. On the one hand it acts as the 
factorization scale, while on the other hand it controls the angular ordering of the partons 
emitted in the evolution |28|, ^9], 0. The factorization scale \i separates the partons 
associated with the emission off both the beam and target protons (in pp collisions) 
and off the hard subprocess. For example, it separates emissions off the beam (with 
polar angle 9 < 90° in c.s.m.) from those off the target (with 9 > 90° in c.s.m.), and 
from the intermediate partons from the hard subprocess. This separation was proved 
in |28], |29|, |30[ and originates from the destructive interference of the different emission 
amplitudes (Feynman diagrams) in the angular boundary regions. 

If the longitudinal momentum fraction is fixed by the hard subprocess, then the limits 
of the angles can be expresseed in terms of the factorization scale /i which corresponds 
to the upper limit of the allowed value of the s-channel parton kr- 

The expression (1) with the survival probability (2) provides the positivity of the 
unintegrated probability f g (x,qT, fJ>) in the whole interval < x < 1. 

Here we deal with the matrix element accounting for the gluon virtualities and 

polarizations. Since it is much more complicated than in the PM we consider only the 

LO of the subprocess gg — > QQ which gives the main contribution to the heavy quark 

production cross section at small x, see the diagrams a, b and c in Fig. 1. The lower 

and upper ladder blobs present the unintegrating gluon distributions f g (y,qiT, ^) and 

fg(x,qzr,IJ>)- 

lr This property is hidden in the conventional parton distributions as qT is integrated up to the scale 
[1. 



The differential cross section of heavy quark hadroproduction has the following form:^ 
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Here s = 2pxPB > ?i,2T are the gluon transverse momenta and ?/* 2 are the heavy quark 
rapidities in the hadron-hadron c.m.s. frame, 



(4) 



|Mqq| 2 is the square of the matrix element for the heavy quark pair hadroproduction. 
It is calculated in the Born approximation of QCD without standard simplifications of 
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the parton model. Contrary to what is mentioned in |L9], the transformation Jacobian 
from x, y to yl, y\ is included in our matrix element. 

We take 

m c = 1.4 GeV, m b = 4.6 GeV , (5) 

for the values of short- distance perturbative quark masses [ 52 , |3"3" |. 



3 Beauty production at Te vat ron- collider 

Eq. (3) enables us to calculate straightforwardly all one-particle distributions as well as 
correlations between two produced heavy quarks using, say, the VEGAS code [Q . 

However there exists a principle problem coming from the infrared region. Since 
the functions f g (y, qir, n) and f g (x, q 2 r, aO are unknown at small values of q\ and qf, 
i.e. in nonperturbative domain we calculate separately the contributions from q\ < Qq, 
Q 2 2 < Ql q{ > Ql and q\ > Q 2 fpj, ||, ||. 

The first contribution (q 2 < Qq, q 2 < Ql) with the matrix element averaged over 
directions of the two-dimensional vectors qy? and q 2 r is exactly the same as the con- 
ventional LO parton model expression. It is multiplied by the 'survival' probability 
T 2 (Ql, fi 2 ). (We assume Ql = 1 GeV 2 in the following numerical calculations.) In this 
contribution the sum of the produced heavy quark transverse momenta is taken to be 
exactly zero. 

The next three contributions (from the domains where one, or both q 2 > Ql) contain 
the corrections to the parton model matrix element due to gluon polarizations, virtuali- 
ties and transverse momenta. Their relative contribution strongly depends on the initial 
energy. If it is not high enough, the first term dominates, and all results are similar 



2 We put the argument of as to be equal to gluon virtuality, which is very close to the BLM 
scheme [pl[; (see also p5[). 



to the conventional LO parton model predictions [[EJ . In the case of very high energy 
the opposite situation takes place, the first term can be considered as a small correction 
and our results differ from the conventional ones. So the highest energies are needed to 
observe (and to study) the fc T -factorization effects. 

In the numerical calculations some uncertainties come from the value of cut-off A in 



Eqs. (1), (2). In particular, in |22| we used 

A = 1 - q T /(j, . (6) 



The angular ordering [^8], ^, Q implies that the cut-off 



A = — ^— . (7) 

In this case one gets the non-zero values of f a (x, Qt, /f) even at the large qx > [A- This 
is especially important at the high energies where the essential values of x and z in 
Eqs. (1), (2) are very small. 

From the formal point of view the difference between (6) and (7) is beyond the 
DGLAP LO accuracy. With the same accuracy we can suggest 

vV 2 + Qt ' 
or something like that. 

Another source of uncertainties comes from the fact that the unintegrated gluon 
distributions are not known with the needed accuracy (from the evolution equation and 
the experimental data). 

Let us compare the results of our numerical calculations with the data on beauty 



production at Tevatron-collider. In Fig. 2a we present the data |35[ on bb pair production 
with pt > Pmin identified by their muon decays, as the function of p m in at y/s = 1.8 TeV. 
The curves of different type show the calculated results with different cut-off A in (1) 
and (2), whereas the upper and lower curves of the same type are obtained with different 
QCD scales, /i 2 = s and /i 2 = s/4, where s is the invariant energy of the produced heavy 
quark pair, s = xys. 

One can see that the scale dependence is negligible here. The dependence on the 
cut-off A is more important. The expressions (6) (dash-dotted curves) and (7) (dashed 
curves) give too small cross section, whereas the cut-off (8) (dotted curves) overestimates 
the data. We obtain the best agreement with the experiment (solid curves), assuming 
the cut-off smaller than (8) and larger than (7), namely 

A = , „ = • (9) 



'fi z + q T + fJ,q T 
The same results for inclusive 6-quark production at large px extracted from muon 



taggeed jets [36| are presented in Fig. 2b. Here again the data are in reasonable 



agreement with our calculations using Eq. (9). 



The rapidity distributions of the produced beauty [|37| are shown in Fig. 3a, b. They 
decrease more steeper than our predictions, but the numerical difference is not large. 

The correlation between the produced heavy quark and antiquark |TT], [38] which are 
more sensitive to the dynamics of the production mechanism are presented in Fig. 4. 
Here the transverse momentum of one beauty quark, pi?, is taken more than some 
value (6.5 GeV in Fig. 4a and 8.75 GeV in Fig. 4b), and cross sections of second quark 
production with p 2 r > Pmin are presented as the functions of p m in- Again, these data 
can be described using the cut-off (9). 

There are the correlations which are more sensitive to the QCD scale value fi in 
comparison with the cut-off value. Some examples are presented in Fig. 5, where the 
distributions over the pair transverse momenta are plotted with the condition that one 
heavy quark has fixed (and large) transverse momentum. In all cases the shapes of the 
curves calculated with fi 2 = s and fi 2 = s/4 are qualitatively different. 

The results for the total cross sections of heavy quark production at FNAL and LHC 
are presented in the Table. These cross sections are essentially larger in comparison 
with the conventional NLO QCD predictions, see, e.g. [Q. It can be important for the 
planning of the experiments at LHC. 



The total cross sections of charm and beauty production calculated in the 
/cr-factorization approach with cut-off (9) and /i 2 = s 





all rapidities 


\y*i\ < i, lx/21 < ' 


V~s 


cc 


bb 
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14TeV 


20.9 mb 


1.16 mb 


2 mb 


142 fib 


1.8 TeV 


3.15 mb 


107 fib 


423 fib 


20 fib 



The azimuthal correlations of heavy quarks, i.e. the distribution over the opening 
angle <fi between the momentum vectors of the produced quarks in the plane transverse to 
the beam axis, are sensitive to the details of the production mechanism. For example, the 
experimental data at fixed target energies are in contradiction with NLO QCD collinear 
approximation [|D|, fi"I ]. Our predictions for bb production at 1.8 TeV and 14 TeV are 
presented in Fig. 6. The same predictions obtained for the large transverse momenta of 
both quarks, pi^r >6 GeV, are also shown. 



4 Heavy quark photoproduction 

In the case of high energy heavy quark photoproduction we should consider (in LO 
QCD) two possibilities, direct production in the photon-gluon fusion, 7p — > QQ, Fig. 7, 
and the resolved production via quark-gluon structure of the photon. In the latter case 
the diagrams are similar to those in Fig. 1, where one proton should be replaced by 
photon. 
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The cross section of the direct interaction can be written as 

,2 



da lp a em e Q f 2 f g (x,q T ,fi) 
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x < Kl _ zf + ^ + ^ +l4 (_L_i_ n , (10 ) 

Di = v\t + m Q , D 2 = (qr - Pit) 2 + m 2 Q , (11) 

where a em = 1/137 and eg is the electric charge of the produced heavy quark. 

The resolved contribution is determined by Eq. (3), where one of the structure func- 
tions f g describes the partons distribution in the proton and the other stands for the 
parton distribution in the photon. 

The energy dependences of the total cross sections of cc and bb photoproduction are 
presented in Fig. 8. The predictions for the cross sections of charm photoproduction 
together with the data [f|^, [43] are shown in Fig. 8a. The averaged cross section value at 
low energies is taken from [ff5| . The data on beauty photoproduction citeHl are shown 
in Fig. 8b. The cut-off value (9) is again in reasonable agreement with the data. The 
values (6) and (7) give too small cross section, as well as NLO QCD approach (44 



These results are in agreement with [45], where different assumptions were used. 

We can see that the cut-off expression (9) gives a better agreement with all the 
presented data in comparison with the other expressions (6)-(8). The fj, 2 dependence of 
all distributions is rather week except of the special case similar to Fig. 5. So, as a rule, 
we will present below the results obtained with the cut-off (9) and /z 2 = s only. 

The one-particle py-distributions, da/dpx, calculated in the fc^-factorization ap- 
proach for HERA energy, are presented in Fig. 9. Here we give separately the sum 
of direct and resolved contributions, as well as the only resolved ones. The resolved 
contributions decrease with pt more fast, that is connected with the finite phase space 
value. At very high energies the resolved contributions will decrease more slow. 

The predicted rapidity distributions, dcr/dy*, of the produced heavy quarks are 
shown in Fig. 10 for all events and for events with px > 6 GeV. The direct produc- 
tion gives the narrow peak in the photon fragmentation region. This peak moves to the 
central region when only quarks with high transverse momenta are registrated, that is 
in agreement with the data [[SJ. The resolved photon contribution dominates in the 
target fragmentation and partly in the central regions. 

The distributions of pair transverse momenta with the condition that one heavy 
quark has fixed (and large in comparison with its mass) transverse momentum are given 
in Fig. 11. Here again, as in Fig. 5, the shapes of the curves calculated with /i 2 = s and 
yU 2 = s/4 are qualitatively different that allows one to determine the scale value. 

The predictions for the azimuthal correlations in heavy quark photoproduction at 
HERA energy are presented in Fig. 12. Their dependences on <fi and on the restriction 
of the minimal transverse momenta of the produced quarks are qualitatively similar to 
the case of hadroproduction shown in Fig. 6. 



5 Conclusion 

We have compared the /c^-factorization approach for heavy quark hadro- and photopro- 
duction at collider energies with the existing experimental data. The agreement is 
reasonable when we take the cut-off (9) in Eqs. (1) and (2) while it does not practically 
depends on the value of QCD scale // in these equations. We present also predictions 
which can be used as an additional test of our approach. Let us note that the different 
expressions for the cut-off A result mainly into normalization factor only, so our predic- 
tions from p2] , |23[| obtained with the cut-off (6) can be used after renormalization with 
the cut-off (9) too. 

Another example of very successful agreement between experimental data and kr- 



factorization approach can be found in Ref. [11] based on different assumptions 



We hope that future experiments will allow to distinguish between various ap- 
proaches. It is necessary to note that the theoretical results concern the heavy quarks 
rather than experimentally observed hadron production. The hadronization leads to 
several important effects, however their description needs additional phenomenological 
assumptions, see e.g. ||7|, |3| ||| |5(], |51[ . 
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Figure Captions 

Fig. 1. Low order QCD diagrams for heavy quark production in pp (pp) collisions 
via gluon-gluon fusion. 

Fig. 2. The cross sections for beauty production with |j/x | < 1, in PP collisions at 
1.8 TeV and their description by the /^-factorization approach with unintegrated gluon 
distribution f g (x, qr, A-0 given by Eq. (1). The values of A in Eq. (1), (2) are taken equal 
to /j,/ ' J(q^ + /U 2 ) (dotted curves), to /i/ ' J(q 2 + fi 2 + qtfi) (solid curves), to aV(<?t + A 4 ) 
(dashed curves) and to \—q T /ji (dash-dotted curves). The upper curves in the left-hand 
side of Fig. 2a correspond to the scale value in Eqs. (1) and (2) /i 2 = s/4, and the lower 
curves to the scale value fi 2 = s. In the case of Fig. 2b upper curves correspond to 
fi 2 = s and lower ones to /z 2 = s/4. 

Fig. 3 The rapidity dependences of beauty production detected by decay muons with 
Pt > 5 GeV (a) and p^ > 8 GeV (b) at 1.8 TeV. The curve types are the same as in 
Fig. 2. 

Fig. 4. Semi differential bb cross sections at \yi\ < 1, |j/2 j < 1, PiTmin = 6.5 GeV (a) 
and 8.75 GeV (b) and their description by the //^-factorization approach with uninte- 
grated gluon distribution f g (x, qx, /i) given by Eq. (1). The curve types are the same as 
in Fig. 2. 

Fig. 5. The distributions of the total transverse momentum p pa ir for bb production 
in pp collisions at ^/s = 1.8 TeV (a) and 14 TeV (b) and with the transverse momentum 
of one heavy quark equal to 20 GeV, calculated with the unintegrated gluon distribution 
f g (x,qT, A-0 given by Eq. (1). The curve types are the same as in Fig. 2. The curves 1 
are calculated with scale /i 2 equal to s, and the curves 2 with the scale /i 2 = s/4. 

Fig. 6. The azimuthal angle distributions of b quarks with |y*| < 1 and jy^! < 1 
produced in pp (pp) collisions at 1.8 TeV (solid and dotted curves) and at 14 TeV 
(dashed and dash-dotted curves) for all px of the produced quarks (solid and dashed 
curves) and for both Pit,P2T >6 GeV (dotted and dash-dotted curves). 

Fig. 7. Low order QCD diagrams for heavy quark production in direct ^p interactions 
via photon-gluon fusion. 

Fig. 8. Total cross section of heavy quark photoproduction. The A values in Eq. (1), 
(2) are taken equal to aVv (<?t + A* 2 ) (dotted curves), to ^/ J(q 2 + fi 2 + qt^) (solid curves) 
and to [x/(qT + /i) (dashed curves). The curves 1 are calculated with scale fi 2 equal to 
s, and the curves 2 with scale fi 2 = s/4. 

Fig. 9. The calculated pr-distributions of charm and beauty photoproduction at W 
= 200 GeV (solid and dashed curves) with A in Eq. (1), (2) equal to fi/J (?t + /i 2 + Q'tA*) 
and scale fi 2 = s. The resolved photon contributions are shown by dotted and dash- 
dotted curves, respectively. 

Fig. 10. The calculated rapidity distributions of charm (a) and beauty (b) photopro- 
duction at W = 200 GeV for all px (solid curve) and px > 6 GeV (dashed curve) with 
A in Eq. (1), (2) equal to /!/ ' J{qr + /i 2 + 9T/-0 and scale a* 2 = s. The resolved photon 
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contributions for the presented two cases are shown by dotted and dash-dotted curves, 
respectively. 

Fig. 11. The distributions of the total transverse momentum p pa i r for cc and bb 
produced in 7p collisions at W = 200 GeV, and with the transverse momentum of a 
one heavy quark equal to 10 GeV , calculated with the unintegrated gluon distribution 
f g (x, qr, /i) given by Eq. (1). The values of A in Eq. (1), (2) are taken equal to /i/y (<?t + 
/-t 2 + 9T/-0- The curves 1 are calculated with scale /i 2 equal to s, and the curves 2 with 
the scale /z 2 = s/4. 

Fig. 12. The azimuthal angle distributions of c and b quarks with all rapidities 
produced in 7p collisions at 200 GeV for all px of produced quarks (solid and dashed 
curves, respectively), for P\t,P2T >2 GeV for charm production (dotted curve) and for 
Pit,P2T >5 GeV for beauty production (dash-dotted curves). 
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